Preparation of atomic velocities by bound-state to resonance conversion 
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A procedure is proposed to control the average and width of the velocity distribution of ultra-cold 
atoms. The atoms are set initially in a bound state of an optical trap formed by an inner red 
detuned laser and an outer blue detuned laser. The bound state is later converted into a resonance 
by a suitable change of the laser intensities. An optimal time dependence of the switching process, 
between the sudden and adiabatic limits, adjusts the final translational energies to the Lorentzian 
shape of the resonance state. 
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I. INTRODUCTION 

With the advent of laser cooling techniques, the tradi- 
tional velocity selection or preparation methods have 
to be substituted, due to the increasing importance of 
gravity and the quantum nature of translational motion, 
by new methods based on Doppler sensitive stimulated 
Raman transitions^, Q, coherent population trapping 
into a dark state or Bragg diffraction by a moving, 
periodic, optical potential IE1E0IEI3- To complement 
these methods and overcome some of their limitations 
, it is worthwhile to explore other approaches based on 
different physical mechanisms. 

Fabry-Perot (FP) matter-wave interferometers realized 
with detuned lasers or microwave cavities have been also 
proposed to provide coherent atomic velocity selection 
or trapping 0, 0, 0, 0, Moreover, the trans- 

mission behavior of a Bose-Einstein condensate [isl IT^ 
through a double barrier in a waveguide has been de- 
scribed in QJl, and through an optical lattice in In 
a recent paper we have explored the fundamental 
limits of a matter-wave Fabry-Perot optical device made 
of two blue-detuned laser barriers and a red-detuned laser 
well, for selecting both the average and the width of the 
atomic velocity distribution. The basic control knob was 
the well-depth, which lets modify the resonance energy. 
It was theoretically and numerically demonstrated that 
this method may produce arbitrarily small velocities but, 
since it is based on filtering the incident velocity distri- 
bution with a resonance peak of the transmission proba- 
bility, the resulting fraction of transmitted atoms may be 
very low and will depend strongly on the incident state. 

The present work describes a modified approach aimed 
at a more efficient preparation and control of the average 
and width of the final velocity distribution. In common 
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with the transmission resonance method, it may also be 
implemented with red and blue detuned lasers, a red de- 
tuned laser forming an inner well and a blue detuned laser 
for an outer barrier. The working principle is now that 
the atom is initialy prepared in a bound state. Then, 
the bound state is converted into a resonance by a suit- 
able change of the laser intensities, so the atom will leak 
out and move asymptotically with the desired velocity 
distribution. This procedure will be descriped in detail 
in Section ^ A necessary requirement for our method 
to be of any practical use is the knowledge of the de- 
pendence between laser parameters and the resonances. 
This relation will be examined in Sec. IIIII An optimal 
time dependence of the switching process, between the 
sudden and adiabatic limits, will de described in Section 
IIVI which adjusts the asymptotic kinetic energies to the 
Lorentzian shape of the resonance state. Final conclu- 
sions and comments will be provided in Section Ivl 



II. BOUND-STATE TO RESONANCE 
CONVERSION 

For simplicity, we shall assume a one dimensional (ID) 
model corresponding to the effective ID atomic motion 
in a narrow waveguide, and a "square" shape for each 
laser intensity, see Fig. ^ although similar results can be 
achieved for smoother profiles The infinite wall at 
the origin is also a simplifying feature of the model, but 
it is not strictly necessary. In particular, one could also 
use two finite barriers, one at each side of the well ^Jl; to 
represent the radial potential profile of a cylindrical con- 
finement with free atomic motion or weak confinement 
in the axial direction. At this stage we also assume a 
simple single-atom or independent-atoms framework de- 
scribed by the Schrodinger equation and disregard non- 
linear effects that could be incorporated within a mean- 
field treatment as in 11]. 

The starting point of the velocity preparation process 
is a laser configuration which holds only one bound state. 
It is assumed that the atom can be prepared in this 
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ground state (see Fig. ^a)). Several possibilities ex- 
ist to prepare that initial state: for example, the original 
trap could hold more than one bound state; in that case 
an arbitrary trapped atomic state overlaps with several 
of them, but the trap may be modified to hold one bound 
state only so that the wave component in the continuum 
subspace is eliminated by its evolution away from the in- 
teraction region. More sophisticated and efficient meth- 
ods without losing atoms may be based on ground-state 
cooling using resolved-sideband transitions Push- 
ing up the potential well later on, the ground state will 
eventually become the only bound state, thus realizing 
our starting point objective. 

Once the initial state of Fig. ^a) is formed, the po- 
tential well is moved upwards by decreasing the intensity 
of the red-detuned laser, i.e. is decreased. In ad- 
dition, the intensity of the blue-detuned laser can also 
be changed. This is represented in Fig. 1(b), where 
the potential switch has been performed suddenly with 
respect to other relevant time scales. The consequence 
is that the bound state becomes, for a final well depth 
shallower than a threshold value, a "resonant state" . As 
it is well known, resonances may be regarded as quasi- 
bound states associated with poles of the S'-matrix in the 
lower half-momentum plane; they can be linked continu- 
ously with bound states (poles on the positive imaginary 
axis) by varying the potential parameters. An impor- 
tant difference though, is that bound states are in Hilbert 
space and normalizable, while Gamow (resonant) states 
are not, since they increase exponentially at large dis- 
tances from the potential center. The normalized state 
achieved by shifting the well bottom, as in Fig. Q^b), is 
thus not a true Gamow state, but it will share approx- 
imately some of its properties, in particular its decay 
rate, the basic Lorentzian shape in energy space and its 
coordinate-space form in the potential region. We insist 
that this agreement is necessarily a partial one. 

After the switching process, the atom will leak out 
(see Fig. QIc)) having a given (total) energy distribution. 
Note that the energy distribution calculated at the end 
of the switching process, i.e., at a time when the atom is 
still interacting with the trap, is equal to the kinetic en- 
ergy distribution of the released atoms at asymtotically 
large time, as it follows from energy conservation. There- 
fore at a sufficient large time, the atom will move with 
the desired velocity distribution. 

Two limits consisting on sudden or infinitely slow well 
switching may be considered: (a) A sudden well shift pro- 
duces a state with contributions from higher resonances. 
They will lead to perturbations with respect to the ideal 
velocity distribution which will affect the short time de- 
cay behavior; (b) The opposite limit of very slow switch- 
ing implies a different problem: since the pole motion in 
the complex momentum plane up to the final desired res- 
onance position is slow, a continuum of other resonances 
are excited. They will have a decay time larger than the 
one desired, thus inducing a deviation with respect to the 
desired exponential decay rate, in this case due to a bias 
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FIG. 1: Schematic representation of the velocity preparation 
method. The left barrier represents an "infinite wall". 



towards slow components. We will see in Section Hvl that 
an adjustment of the switching time may avoid the per- 
turbations of the fast and the slow processes and produce 
an excellent agreement with the Lorentzian shape of the 
energy distribution. 

An implementation of the proposed velocity prepara- 
tion method will require the knowledge of the dependence 
between well/barrier parameters and the resonances such 
that the atom leaving the trap will have the desired ve- 
locity distribution. This dependence can be achieved ex- 
perimentally or theoretically, as it will be described in 
the next section. 



III. CONFIGURATIONS AND 
CORRESPONDING RESONANCES 

We consider the model based on a well and a barrier 
with variable depth/height represented in Fig. ^ The 
effective potential is assumed to take the initial and final 
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In the final configuration, see Fig. ^b), the stationary 
states of a single ultra-cold atom moving along the x 
direction will satisfy 



2m dx^ 



where Ek = For the calculations we have used the 

mass of ^^Na. The scattering states will have the form 
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where q = yk^ + 2mVw/h and q' = yk^ — 2mVb/h 

(for the rest of this section we omit the superscript ) . 
q and q' have a branch cut in the p-plane joining the 

two branch points at ±i^2mVw/h^ and ±\J 2mVh/h? , 
respectively. The different coefficients are obtained from 
the matching conditions at x = 0, a; = d, and x = d + b. 
The resonances and bound states can be calculated from 
the poles of the S'-matrix in the complex fc-plane. They 
are solutions of the equation 

n{k) -{k^q')[q + e^^'''i{q-q') + q'] 

+ e^'^'t' (k + q') [q-q' + e^'^'^iq + q')] = 0. 

The corresponding roots in the upper half-imaginary axis 
are the bound states of the system, while the roots in the 
fourth and third quadrant are resonances and antircso- 
nances, respectively. 

An alternative way to find the resonances is to look for 
jumps of the phase shift S{k), 



log[5(fc)]. 



or the peaks of the Wigner delay time At, 



At{k) = 2h 



d6{Ek) _ 2md6{k) 



dEk 



hk dk 



(1) 



This may be easier than determining the poles by ana- 
lytical continuation of S{k). In the Breit- Wigner regime 
of isolated and sharp resonances, the information about 
the resonances given by these peaks may be straightfor- 
wardly related to the poles in the complex fc-plane. If 
Erea = Efi — iT /2 is the complex energy obtained from 
the phase shift (central position En and width at half 
height r), then 
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FIG. 2: Barrier height Vb versus well depth Vw (upper panel) 
and, energy width F versus Vw (lower panel) for constant 
resonance energies Er: Er = 53.391?i/s (solid lines) and 
Er = 7.422?i/s (dashed lines). 6 = 10 /im, d = 5 /im. 



where k^es = fci + ifc2 is the corresponding complex wave- 
number for which ^{kres) — 0. 

Different combinations of well depth Vw and barrier 
height Vb lead to the same resonant energy Er. In Fig. 
|2ta) we have plotted the curves for two different resonant 
energies Er with 5 = 10 /xm, and c? = 5 ^m. Along a 
curve for a fixed Er, the resonance width will change. 
This change is plotted in Fig. EJb). Notice that for small 
values of the well depth V^, the barrier height Vb that 
keeps the resonant energy Er fixed is almost constant. 
In contrast, F changes drastically for small changes of V^ 
in shallow wells. 



IV. RESULTS FOR DIFFERENT SWITCHING 
TIMES 

Let 'fo{x) be the ground state in the initial potential 
configuration of Fig. H^a) with 1/^"'* = 350a/s, V/"** ^ 
4007i/s. For the destination resonance in Fig. Q^b), we 
choose the lowest one corresponding to V^*" = 100?i/s 
and V"j/^*" = 200?i/s. In this case, the resonant complex 
energy is Eres = (134.509 — il.2n)h/ s. To characterize 
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FIG. 3: Delay time (solid line, left axis), versus incident 
energy of the chosen resonance (see text); the thick dashed 
line corresponds to a fitting with a Lorentzian profile. Energy 
distribution P{E) (circles, right axis) of the state which is 
initially the ground state for 1/^"^* = 350?i/s, V/"'' = 400;i/s, 
in the final configuration T<f'" - 
d = 5 fim, fo = 10 ^im. 



= lOOft/s, V^""' = 200;i/s, 



this resonance we have plotted in Fig. O the delay time 
At, see Eq. versus the incident energy. Also shown 
is the Lorentzian fitting. 

If we move suddenly the bottom of the well making it 
shallower until the initial bound state overlaps strongly 
with the desired resonance (Fig. 1(b)), the wave- function 
will evolve in time, and the atoms will leak out through 
the barrier (Fig. 1(c)). We assume that the ensemble 
of non-interacting atoms satisfy during this process the 
one-dimensional time dependent Schrodinger equation. 



ih 



dtjj{x, t) 
'dt 



2m dx^ 



iP{x,t), (2) 



where ij^i^j 0) = tpo{x). As described above, our main ob- 
jective is to achieve a distribution as close as possible to 
the Lorentzian distribution associated with a Gamow res- 
onance of the final potential configuration (Figs. 1(b) and 
1(c)). Fig. 13 shows also the resulting total energy dis- 
tribution P{E) = y/ I (-^fc^ g) lyo) P of the wavepacket 
in the new potential configuration. {P{E) coincides with 
the kinetic energy distribution when the packet moves 
away from the potential region. 25] ) After such a sudden 
process, several resonances are excited as it may me seen 
in different ways: note in particular that the Lorentzian 
of Fig. 121 is not normalized. This means that part of the 
norm is in higher resonances. A consequence is the fast 
decay of the non-escape probability 



Pw{t) 



dx \'ip{x,t)\'^' 



at short times in Fig. 0] (solid line). In other words, 
with the sudden switching a significant fraction of atoms 
is released at early times with "too much" energy. Of 
course, if we discard the early, fast atoms, the decay oc- 
curs finally with the desired rate and energy distribution. 
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FIG. 4: Decay of the probability Pw{t) to find the atom in 
the well versus time for different values of T, see Eq. ©; the 
initial and final potential configuration are given in Fig. |3 



see Fig. 01 However, we may try to produce an ensemble 
without undesired high velocity components. This can be 
achieved by a progressive, rather than abrupt, switching 
of the potentials. 

Let us assume that the potential profile changes in time 
according to the smooth function 

V{t,x) = [^^"(2:)-y^""(x)](l-e-*/^)-f y^™*(2;).(3) 

The sudden change corresponds to T = and the in- 
finitely slow change to T = oo. In Fig. 01 we show the 
decay of the non-escape probability Pw (t) ■ The lifetime 
of the first resonance (calculated from the pole of the S 
matrix) is r = 0.411 s, in perfect agreement with the 
fitting to the exponential decay that dominates after the 
early transients, independent of T. This occurs because 
the final potential configuration is common to all cases 
so that the resonance with the longest life time is the 
same in all cases. Nevertheless, the first transient regime 
varies substantially with T, and for T w r the initial 
decay is slowed down considerably. The best fit to the 
purely exponential decay is found for T « 0.13r. 

Our interest is in the asymptotic and stationary energy 
distribution at large time, too ^ T. In Fig. [3 we have 
plotted the final energy distribution for different transi- 
tion times T. For T = the distribution at the main 
resonance can be approximated by the normalized Breit- 
Wigner Lorentzian form corresponding to the pole of the 
selected resonance, however its height is reduced because 
of the excitation of higher resonances. Increasing T, an 
optimal value is found so that the energy distribution fits 
even in magnitude to the ideal Lorentzian shape. For the 
case studied in Fig. [3 the best fit to the Lorentzian of 
the selected resonance corresponds to T « 0.058t. Note 
that this optimal value of T is different from the one 
that provides the best fit to the purely exponential decay 
(T « 0.13r), compare Figs. ^andO As T is increased 
further, the distribution is deformed, the symmetry is 
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FIG. 5: Final energy distribution for different transition 
times T. The dots correspond to the normalized Lorentzian 
distribution with the parameters extracted from the S-matrix. 
The initial and final configuration are the same than in Fig. 



lost and a distortion favoring lower energies is observed. 



V. CONCLUSIONS AND DISCUSSION 

We have proposed a method to prepare states with well 
defined average velocity and width, based on the conver- 
sion between a bound state and a resonance by changing 
the trapping potential. In atom optics the potentials may 
be realized with detuned lasers. An optimal switching 
time exists so that the kinetic energy distribution of the 
leaking atoms fits the Lorentzian form of the resonance. 

Laser fluctuations may tend to broaden or blur quan- 
tum resonances. Nevertheless, stabilized lasers provide 
effective constant intensities in the time scale of r and T 
(~ 0.01 - 1 s) liO, m 113 so, that our analysis would 
apply to the effective, time-averaged potentials. Further 
theoretical and experimental work is required to deter- 
mine the feasibility of using scattering resonances in op- 
tical traps. 

Finally, the method may be applied to electrons in 
semiconductor heterostructures, where the well-depth is 
modified by potential gate voltages |3 . 
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